PROCEEDINGS 

™ SCIENCE 



Thermal dilepton rates from quenched lattice QCD 



o 

O 



i3 



> 

o 



X 
S3 



O. Kaczmarek, E. Laermann, M. Miiller* 

Fakultat fur Physik, Universitat Bielefeld, D-33615 Bielefeld, Germany 
E-mail: okacz, edwin, mmueller@physik.uni-bielefeld.de 

F. Karsch 

Fakultat filr Physik, Universitat Bielefeld, D-33615 Bielefeld, Germany 
Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
E-mail: karschgquark . phy . bnl . gov 

H.-T. Ding, S. Mukherjee 

Physics Department, Brookhaven National Laboratory, Upton, NY 1 1973, USA 
E-mail: htding, swagato@quark.phy.bnl.gov 

A. Francis 

Institut filr Kernphysik, Johannes Gutenberg-Universitdt Mainz, D-5 5099 Mainz, Germany 
E-mail: f rancis@kph . uni-mainz . de 

W. Soeidner 

Fakultat fur Physik, Universitat Regensburg, D-93053 Regensburg, Germany 
E-mail: wolf gang . soeldner@physik . uni-regensburg . de 

We present new lattice results on the continuum extrapolation of the vector current correlation 
function. Lattice calculations have been carried out in the deconfined phase at a temperature of 
T = l.lTc, extending our previous results at T = 1.45rc, utilizing quenched non-perturbatively 
clover-improved Wilson fermions and light quark masses. A systematic analysis on multiple 
lattice spacings allows to perform the continuum limit of the correlation function and to extract 
spectral properties in the continuum limit. 

Our current analysis suggests the results for the electrical conductivity are proportional to the 
temperature and the thermal dilepton rates in the quark gluon plasma are comparable for both 
temperatures. Preliminary results of the continuum extrapolated correlation function at finite 
momenta, which relates to thermal photon rates, are also presented. 
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1. Introduction 

Current Heavy-Ion-experiments at LHC and RHIC reach energies in a region where the transi- 
tion from normal nuclear matter to a quark gluon plasma occurs. In this phase thermally produced 
dileptons and photons are important experimental observables. To understand their yields, hydro- 
dynamic models for the evolution of the medium need detailed knowledge of the dilepton and 
photon rates as well as the transport coefficients during the evolution of the system [1]. These 
require non-perturbative ab initio lattice QCD calculations as input. 

On the lattice, hadronic coiTclation functions for different particle channels can be calculated at 
zero and at fixed finite temperatures. From the vector spectral function the electrical conductivity 
and the dilepton rates of the medium can be extracted. The spectral function has recently been 
extracted by employing phenomenologically inspired ansatze [2, 3]. 

Our previous work on dilepton rates provided the first results for a continuum extrapolation 
of the vector meson correlation at a temperature of T /T^ = 1.45 and light quark masses [4]. An 
ansatz for the spectral function could successfully be used to describe this dataset, thereby the 
dilepton rate could be obtained at a first relevant temperature. Together with results for the electrical 
conductivity at different temperatures 1.16 < T /T^ < 2.98 [5] this motivated the systematic study 
of the temperature dependence of the vector spectral function on the lattice. For results using two 
dynamical flavors at finite lattice spacings see [6]. 

The work presented here extends the continuum extrapolated quenched calculations to a sec- 
ond relevant temperature T /T^ = l.l. 

2. Thermal vector correlator and spectral function 

2.1 Vector correlation function 

The Euclidean time two-point correlation function G(t, p) of the vector current is a quantity 
directly accessible in lattice QCD calculations. 



Only contributions of quark line connected diagrams ai^e included. Disconnected diagrams 
cause a high numerical effort, but aie expected to be small in the high temperature phase of QCD 
[7, 8]. The correlation function directly relates to the spectral function via 



Here pu denotes a sum over the spatial components, while poo denotes the time-like compo- 
nents. The full vector spectral function is denoted by py = poo + Pu- As lattice QCD calculations 
are carried out at a fixed temperature, the explicit temperature dependence of p is dropped and 
observables are usually given in units of the temperature T. 

2.2 Ansatz for the spectral function 

The time-like component of the vector coiTclator Goo and thereby the corresponding transfor- 
mation of the spectral function poo is related to the quark number susceptibility Xq- Since the quark 
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(2.2) 
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number is conserved, the correlator is constant in (here EucUdean) time, Goo(Tr) = —XqT and its 
spectral representation is given by a delta function 

poo(ft>) = -2tiXc,0}5{(0). (2.3) 

The spatial components of the spectral functions increase quadratically for large values of ft), 
in the free field limit for massless quarks to 

pj'^^ico) = 27iT^(o5((o) + ^(0^-tanh((o/4T). (2.4) 

271 

In this limit, the delta peak in the spatial and the time-like component of the spectral function 
cancel. In the interacting theory however, this is not the case: The time-like component maintains 
a delta peak since it is linked to the conserved current. In the spatial component however the delta 
peak is smeared out and expected to be described by a Breit-Wigner peak [9, 10, 1 1, 12] 

pr"^-(«) = X,CBW-^^^ + (1 + K)^cohanHco/4T). (2.5) 

This ansatz leaves three parameters, the strength (cbw) and width (F) of the Breit-Wigner peak 
as well as k, which accounts for the deviation from free theory. The relation of this ansatz to the 
correlator obtained on the lattice is given by (2.2). 

The fits are not performed directly to the correlation function G„ but to a set of two ratios: 
The correlation function is normalized by the quark number susceptibility (as given in (2.3)), re- 
sulting in a dimensionless quantity independent of renormalization constants. It is also normalized 
by the free field correlation function Gy^^{zT), yielding a smooth function that does not fall off 
over multiple decades like the correlation function. Furthermore, due to asymptotic freedom, the 
con^elation function approaches the non-interacting limit at asymptotically small distances. The 
spectral function is thereby fitted to reproduce 

G.(Tr)/Goo ^2.6) 



Having obtained the spectral function, relevant properties of the medium can be calculated, 
e.g. the electrical conductivity as 

a ^ a(r)/Ce.„ = 2;^,CBw/(3F) (2.7) 

1 6 £0^0 (01 

where Cem is given by the elementary charges Q of the quark flavor / as Cgm = ^^'^ 
thermal production rate of dilepton pairs as 

dW 5a^ \ , , 

Pu{(0,p,T). (2.8) 



dcod^p 547r3 ftj2(f.ffl/r_i 

2.3 Thermal moments 

Using the ansatz above to perform a three-parameter fit of the spectral function, the Breit- 
Wigner parameters are most sensitive to the low CO region, thereby to the large distance behavior of 
the coiTclator. 
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No 




K 


l/a[GeV] 


a[fm] 


#conf 


32 


96 


7.192 


0.13440 


9.65 


0.020 


314 


48 


144 


7.544 


0.13383 


14.21 


0.015 


315 


64 


192 


7.793 


0.13345 


19.30 


0.010 


242 



Table 1: Summary of simulation parameters 



This region of the correlator can be further constrained by calculating the curvature of the 
correlator for large Euclidean distances xT . The thermal moments 



(„) 1 d'^GniTT) 
(jf ~ 



1 Z"" dco /- (0\" Ph((0) 
Tr=i/2 "'"'0 271 Vr/ sinh(G)/2r) 



" n\ d{xTY 

are given as the Taylor coefficients of the correlation function expanded around the midpoint 

GH{TT) = Y^G^^''\--xTf\ (2.10) 

For the infinite temperature, free field limit, the thermal moments can be calculated analyti- 
cally. For the analysis, the ratios of interacting to free midpoint subtracted correlation functions, 



have been calculated, where /?y = Gy ^ / G^T^ ■ The first two terms (up to quadratic) are then 
used in the fit. 

For the 1.45rc dataset, restricting the spectral function fit to also reproduce these first two 
thermal moments allowed to obtain a stable fit. Since the new lATc dataset allows for a more 
precise continuum extrapolation down to a distance xT^^ = 0. 15, which is smaller than the distance 
'^Tmin = 0.25 reachable at 1.45rc, we find that stable fits at this temperature (giving results within 
the systematic errors of this analysis) can also be performed without including the thermal moments 
as an additional constraint. 



3. Results for correlators and spectral function fits 

Vector correlators are obtained on lattices of size A^^ x Nt, where for the new dataset at 1.1 
three fixed lattice spacings {aT)^^ =Nt = 32,48,64 where chosen for the continuum extrapolation. 
For each size jS was set to match the fixed temperature and K was tuned to light quark masses (see 
table 1). From the previous study at 1.457^., finite volume effects are found to be small for a fixed 
aspect ratio Na/Nx = 3, so the lattice sizes were set to 96^ x 32 , 144^ x 48 and 192^ x 64. 

Cut-off effects are removed by a continuum extrapolation of the coiTclators: Discretization 
eiTors of non-perturbatively clover-improved Wilson fermions have a quadratic eiTor in the lattice 
spacing allowing to extrapolate the correlators in a^, corresponding to l/N^ at fixed temperature 
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Figure 1: Left: Continuum extrapolation of the correlator ratio Gv{TT)lGy^^'^'^'{'iT) and 
Gv{tT) / Gy^^'''°"^ {xT) at different time-like separations xT . Cutoff effects are under control for 0.125 < 
XT < 0.5. Right: Correlator ratio Gv(Tr)/G{).''''™"' (Tr) for the three lattice spacings and continuum ex- 
trapolation. 



T . For these extrapolations, the ratios (2.6) of the interacting correlators to free theory are used. 
On the coarser lattices these ratios are spline-interpolated to provide every spacing xT present 
in the finest correlator. As can be seen in figure 1, calculating ratios to the free continuum as 
well as to the free lattice correlators give the same continuum extrapolation for Euclidean time 
separations xT > 0.125 and the results show that lattice cutoff effects are under control in the 
region 0.15 < tJ < 0.5. 

3.1 Fitting the spectral function 

After the continuum extrapolation of the vector correlation function has been obtained, the 
first and second thermal moments are computed by fitting the curvature of the midpoint subtracted 
correlator (2.11). Both the correlator ratio Gy {xT) / XqGy'^'^ {xT) and the ratio of the first thermal 
moments serve as input and constraints for a fit of the parameters cbw,^, ^ in the spectral function 
ansatz (2.5) to reproduce the correlator (2.2). As can be seen in figure 2 (right), the fitted spectral 
function reproduces the correlation functions. The thermal moments constrain the small CO region, 
so errors of the spectral function fit descrease with increasing xT. 

3.2 Breit-Wigner with truncated continuum 

As in the analysis of the 1.457^ dataset [4], to study the systematic uncertainties of our Ansatz 
for the spectral function the low energy structure of the spectral function was studied by smoothly 
cutting the continuum contribution at a frequency (Oq with a width of A(b 



pjr ico) = X,CBW ^2 4^r/2)2 + + K)^cohanh{w/4T)eiwo,Aa>) (3.1) 



with 0((Oo,A«)= (l+e("o-'"')/^^-] 
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Figure 2: Left: Ratio of midpoint substracted correlators (Gy(Tr) -G[,°V(Gv'''(Tr)-G[,'''*''), including 
a fit of the curvature (2.11) and thereby the thermal moments. Right: Fit of the continuum extrapolated 
correlator to (2.2) with the spectral function ansatz (2.5). 



The fit with (3. 1) is performed for a range of values (Oq and A^. The systematic errors given in 
figure 3 correspond to the minimal and maximal values cbw found in these fits with ;^^/d.o.f < 1.1, 
thereby giving a minimal and maximal electrical conductivity (2.7) within this framework. 

From the corresponding spectral functions the dilepton rates shown in figure 3 (right) are 
calculated (2.8). Within the current systematic uncertainties, the dilepton rates and electrical con- 
ductivities in units of T agree between \.\Tc and 1.457^. This linear temperature dependence of 
the electrical conductivity is in agreement with the results of [13]. 

In order to study systematic errors, also a variation of the ansatz (3.1) motivated by Operator 
Product Expansion [14] can be used. First tests show only small deviations that are well within the 
systematic errors between both version of the ansatz. 



4. Outlook: Correlation functions at finite momenta 



Extracting the spectral function at vanishing momentum p = gives access to the dilepton 
rates and the electric conductivity. For finite momenta p the correlation and spectral functions 
split into a longitudinal and transversal part, where the transversal pr relates to the photon rate: 

dR^^ pj{co = \plT) 

d^p exp(w/r)-l ■ ^ ' ' 

On the lattice, only a finte number of discrete momenta are accessible, the momenta are given 
by the aspect ratio Na/Nx and a vector of integer numbers k as p/T = 2k -k ■ Nt/No- All I.IT^ 
correlator calculations have been performed on the same aspect ratio of A/^cj/A^t = 3, so a continuum 
extrapolation of the coiTclator at fixed momenta can be performed. As can be seen in figure 4, this 
extrapolation is well behaved as in the zero momentum case and gives precise results for the vector 
correlation function at finite momentum p in the continuum limit. Work on the extraction of the 
corresponding spectral functions and the determination of the photon rate is in progress. 
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Figure 3: Left: Spectral function as fitted to the vector correlation function. Two lines are shown for 
each temperature as a result of the systematic error estimation, given by the two pu with the minimal and 
maximal cbw that can be obtaind via the Breit-Wigner + truncated continuum ansatz (3.1), while maintaining 
a x/d.o.f < 1.1 in the fit. Results from the hard thermal loop resummation scheme (HTL) are also included 
in the plot [15]. Right: Thermal dilepton rate calculated from the spectral function (2.8). 
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Figure 4: Continuum extrapolation of the logitudinal (left) and transversal (right) vector correlation function 

-•fre 
'V 



Gv I G^^'^ at non-zero momentum. 



5. Conclusion 

Lattice calculations of the vector correlation function have been performed for two temper- 
atures in the deconfined phase of QCD. Calculations at different lattice spacings were used for a 
successfull continuum extrapolation to remove cut-off effects. Using a phenomenologically moti- 
vated ansatz, the vector spectral function to the continuum-extrapolated correlation function was 
obtained, giving access to the dilepton rate and the electric conductivity of the medium at the 
given temperatures. Within current systematic uncertainties, the electric conductivity divided by 
temperature and the thermal dilepton rates are are compatible at L 1 and 1 .457^. A continuum ex- 
trapolation was also performed for the correlation functions at a range of finite momenta, opening 
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the possibility to study the photon rate in the future. 
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